% Funcion simular
% Parametros de entrada:
%	- x0 son las condiciones iniciales
%	- modelo es el ivp que se desea simular
%	- t0, tfin es el intervalo de tiempo
%	- params son los parametros que recibe el modelo
%	- h es el tamanio del paso de integracion
%
% Parametros de salida:
%	- ths es el angulo
%	- ts son los tiempos asociados a todos los demas vectores
%	- theta_r

function [ths, ts, theta_r] = simular(x0, t0, tfin, modelo, params, h);
	[ts xs theta_r] = rk4(modelo, params, t0, tfin, x0, h);
	ths = xs(:, 1);
end
